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Abstract 

In the theory of superfluidity and superconductivity, a jump of 
the free energy was discovered theoretically and was naturally called 
a zeroth-order phase transition. We present an example of an exactly 
solvable problem in which such a phase transition occurs. 

1 Introduction 

The Nuclear Power Plants (NPPs) cannot be stopped comletely after their 
resources have been exhausted. And this fact of "long-living" means that 
the well-known decay microprocesses develop in the depths of the NPPs 
after their service time is over or after an accident occurs. We must study 
these microprocesses from the standpoint of macrophysics, i.e., we must use 
the averaged thermodynamical approach. 

At present, we have survived several accidents in NPPs the greatest of 
which is the accident in the 4th block of the Chernobyl NPP. The resources 
of the Leningrad NPP will soon be over (and this NPP has already issued a 
radioactive release). What shall we do further with this NPP? It was decided 
to continue its operation time. 

The service times of the French NPPs will also come to the end in the 
near future. How will they behave after the reactor is shut down? It is 
impossible to perform such experiments in laboratories. 

This is a completely new phenomenon, and we face a problem which has 
never been posed before: to understand what the "black box" (an NPP very 
slowly dying after an accident or after the end of its service time) contains. 

In the spring of 1986, I was invited to collect a group of mathematicians 
(and to be the head of this group) whose task was to eliminate the effects of 
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the accident in the 4th block of the Chernobyl NPP. The enormous respon- 
sibility, very short terms, unexpectedness of some new phenomena observed 
for the first time demanded a huge strain of our forces and completely turned 
over my understanding of mathematics as a science. 

In 1987 we published our report in the form of book [T], but we did not 
find answers to the main questions. 

It was clear that the modern mathematics must be revised in order to 
generalize and intensify the thermodynamics, the theory of probabilities, and 
the mathematical physics, i.e., the branches of science required to solve the 
problems arising in NPPs after accidents. 

Since we cannot get inside a broken reactor and perform measurements 
there, we obtain only average results of measurements, i.e., thermodynamical 
results. Inside the broken reactor, there may be no radiation at some places, 
and dangerous spots are spread at random. To explain several factors, we 
have to revise and generalize the thermodynamical science with regard to un- 
expected physical phenomena. Since the processes of nuclear decay continue 
in an NPP after an accident or after its service time is over, we must con- 
sider these processes from the thermodynamical standpoint assuming that 
the "service time" of the NPP is infinitely long. 

In general, it is already somewhat to late to revise the basic notions of 
the classical thermodynamics, because, after the superfluidity was discov- 
ered, the idea that all the molecules stop their motion at the temperature 
zero was refuted by the experiments. Therefore, it is necessary to study 
the superfluidity from the physical standpoint and from the standpoint of 
physical-and-mathematical generalizations. It is necessary to generalize the 
remarkable results obtained by Landau's and Bogolyubov's schools so that it 
were possible to use them in the seemingly opposite situations arising after 
accidents in NPPs. 

The works of physicists of genius intuition must later be generalized from 
the mathematical-and-physical and, perhaps, philosophical standpoint by 
theorists of such a high level as von Neumann, Wigner, and others. 

In his letters, A. N. Kolmogorov cites P. L. Kapitsa's opinion about su- 
perfluidity. Kolmogorov writes: "P. L. (Kapitsa - Victor Maslov) reported 
the following: already before the war, he (Kapitsa - Victor Maslov) per- 
formed a series of experiments for verifying Landau's theory. The results of 
these experiments (which had not been published) absolutely contradicted 
all Landau's constructions. . . . Obviously, the publications were delayed on 
purpose, so that Landau could not correct his theory (because, as is known, 
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'he could prepare any theory justifying any given data'!)" |2j, p. 594. 

In fact, to give an instantaneous response to any new experimental data 
is a talent and a step towards the truth. And to relate the concept of super- 
fluidity, the classical thermodynamics, the theory of averaging (the theory of 
probabilities), and the data of accidents as a whole, this is already a problem 
for a mathematical physicist. 

I had a powerful stimulus to establish this relationship, because the prob- 
lems that I could not solve at the time of the Chernobyl accident have been 
bothering me for the last 20 years, and the emotional experience I went 
through when signing the permission to close the sarcophagus without being 
fully confident have been torturing me even in my sleep. Each of the regular 
releases in the NPP returned me to these problems. 

To solve this problem, we must combine the statistical physics and the 
science of its averaging in physics (the thermodynamics), as well as the math- 
ematical statistics and the science of its averaging (the theory of probabili- 
ties), as a single science. 

But the modern theory of probabilities in not sufficient to deal with the 
processes in NPPs. If the probability of an accident is negligibly small, 
but the effects of the accident are grandiose and disastrous, then even an 
"infinitely" small probability cannot be neglected. 

The second point is the following. Let us consider a mathematical ex- 
ample about the heads and tails play. The probability of a head or a tail 
at a given moment is 1/2 independently of the prehistory of heads and tails 
till the present moment. But, as the famous English physicist R. Peierls 
wrote about the roulette play (black and red instead of heads and tails) , it is 
doubtful that "there would be at least one player who, once taking the risk of 
black, will continue to stake constantly the same number . . . the probability 
of the appearance of such a long "black" series is negligibly small so that it 
must be ignored completely" [3], P- 126. How does this agree with the fact 
that the probability is independent of the prehistory? 

These considerations are very close to the notion of "general position" in 
topology. So, for example, the "general position" of a curve on the plane with 
respect to the projection on the x-axis is the position of the curve such that 
the "cardinality" of the set of nearby curves with the same properties is much 
greater than the "cardinality" of the set of curves with different properties. 
And this construction is similar to the considerations given by Boltzmann in 
answer to the Zermelo objection: if one state is arbitrarily chosen from the 
set of all possible initial states, then it is most incredible that this will be a 
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state with decreasing entropy. This notion, which is close to that of stability, 
does not belong to the mathematical theory of probabilities, but it must 
belong to the new revised "theory of probabilities" because it corresponds to 
the class of problems mentioned above. 

It should be noted that the common mathematical notions such as "nec- 
essary conditions" and "sufficient conditions" are not sufficient for solving 
the specific problems arising in NPPs after accidents. For example, it is pos- 
sible to write necessary conditions for constructing the sarcophagus so that a 
release will necessarily occur, but it is more than difficult to write "sufficient 
conditions" for ensuring that no release will occur. 

How can we choose the most reliable "stable" construction from a set of 
equivalent constructions? This problem is close to the Peierls "law" and to 
the notion of "general position." 

The general thermodynamical approach with the relationship between 
the service time and the time of relaxation to the quasi-equilibrium state, 
(i.e., the asymptotics in general rules of averaging) taken into account shows 
that there is something common between such phenomena as the spouting 
of superfluid helium in Allen- Jones experiments in 1938 and the radioactive 
releases in NPPs. This is the zeroth-order phase transition, which has not 
been noticed by theoretical physicists jU E]- in this phase transition, not 
only the entropy (the intrinsic energy), but also the thermodynamical Gibbs 
potential, has a jump. The instant at which this can occur can be determined 
from the corresponding equations of state. This effect essentially differs from 
the well-known picture of the first-order phase transition (cf. Fig. 4 in [5J). 

A hydrodynamical rather slowly process was also be discovered and stud- 
ied in the rough approximation in pQ. 

It was a common opinion that fuel rods will necessarily melt everything 
under them and go through the ground into the subsoil waters as an iron 
goes through ice. This could cause an ecological disaster. To avoid this, even 
before our calculations were finished, a massive concrete basement was built 
under the debris in order to keep the fuel rods from going under the ground. 
But this enormous work was, in fact, done for nothing, since the calculations 
showed that this undesirable effect could be avoided. 

The author and the late V. P. Myasnikov found that this problem can be 
solved by using the Darcy equation, i.e., the equation of filtration. To solve 
this problem, we performed several experiments and numerous mathematical 
calculations. It turned out that, even in the case of filtration, because of the 
difference in pressure, a small blow from beneath decreases the temperature 
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under the fuel rods. This phenomenon is similar to the chimney effect when 
the chimney grate is not heated, although the temperature of the burning 
coal is very high. This was confirmed in practice, none of the fuel rods came 
to the concrete basement. 

The three-dimensional calculations taking the construction of the sar- 
cophagus and the hydrodynamical process of "convection" in the filtration 
medium into account were performed by the author more precisely in [UJ |S] 
and in 0, pp. 245-325. 

2 General concepts concerning 

the zeroth-order phase transition 

2.1 Main large parameters in thermodynamics 

Thermodynamics studies steady-state processes in which, independently of 
its initial state, the system comes to a state that remains the same in what 
follows. But if such a state still varies, then this is a thermodynamic variation 
only if the process occurs extremely slowly PU|- in other words, we slightly 
vary the state of the system and wait until it returns to equilibrium and 
again is independent of the initial state. 

Thus, a large parameter, the time, is invisibly present in this case i.e., 
we observe the system at huge time intervals. On the other hand, ther- 
modynamics is the limit macroscopic theory obtained from the microscopic 
statistical physics as the number of particles N tends to infinity. Hence there 
are two large parameters, and many things depend on their ratio. In turn, 
the classical statistical physics is the limit of quantum statistics as h — > (h 
is the Planck parameter). Thus, already three large parameters "meet" here: 
time scale, number of particles, and 1/h. 

In thermodynamics, sufficiently large time periods are also often consid- 
ered, and in these periods only a part of the system comes to equilibrium; 
this is the so-called "local equilibrium." For example, when plasma is heated 
by a magnetic field, the ions begin to obey the Maxwell distribution after a 
large observation period, and the entire system comes to the thermodynamic 
equilibrium only after a significantly large time period. 
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2.2 The Gibbs postulate 



Let us consider the Gibbs postulate, which has the following form in quantum 
statistics. 

Suppose that a system is characterized by the Hamiltonian operator H 
in the Hilbert space (in particular, in the Fock space) $. Suppose that the 
operator H has a nonnegative discrete spectrum A , Ai, . . . , \ n . . . . Then the 
free energy depending on the temperature 6 is determined as 



where 5 n is the multiplicity of the eigenvalue A„. 

We note that one more parameter appears in the Gibbs postulate, namely, 



This is one of the most important facts, since liniju-^oo and liniAr^oo do not 
commute! It turns out that one must first pass to the limit as M — * oo and 
then as N — > oo. 

If we speak about local equilibrium, then, as a rule, the Gibbs formula 
deals with a subset {A n /} of the set {A„}. This must be done in the study 
of the above problem concerning the ions distribution in plasma on heating. 
This is also done in glass and in several other physical problems. 

2.3 "Friction" operator 

At temperature 8 = 0, the Gibbs formula gives the value of the lowest eigen- 
value, which, in physics, is called the "ground state." 

This is a law which has not been stated accurately and which was called 
the "energetic efficiency" by N. N. Bogolyubov and others. The lower energy 
is energetically more efficient. For example, if a system is perturbed by an 
operator V(Q — > $) and this operator is sufficiently small, then the transition 
matrix from the state X n to the state A m is determined by the matrix element 
(^ n V^^), where ^ n and \l/ m are the eigenf unctions corresponding to A n 
and A m , respectively (see [11J). 

The square of this matrix element is the probability of transition from 
the state with energy A n to the state with energy A m . But if A n < A m , then 
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this transition is "not energetically efficient" and hence, as if from the view- 
point of statistical physics and thermodynamics, is unrealizable, i.e., from 
the mathematical viewpoint, it must be set to be zero. Then the transi- 
tion matrix is not a self-adjoint matrix with zero entries above the diagonal. 
This means that the friction is taken into account. A pendulum oscillates 
and stops in the end if the friction is taken into account, i.e., comes to the 
"ground state." The "friction" operator was considered in more detail in 
author's paper |T2] . 

2.4 Phase transitions 

The derivative of e{9) with respect to 9 is called the entropy, and when the 
entropy has a jump at a point #0, it is said that a first-order phase transition 
occurs; when the second-order derivative has a jump, it is said that a second- 
order phase transition occurs, etc. 

But, of course, not the function itself has a jump, but the leading term 
of its asymptotics as iV — > oo. In experiments, this "jump" can sometimes 
be actually extended in time, but we have agreed to consider variations in 
sufficiently large periods of time. In these periods, not only thermodynamic, 
but also dynamical processes can occur. We neglect them and consider only 
the time periods in which the system comes to equilibrium with the "energetic 
efficiency" taken into account [T3] . 

The author discovered the zeroth-order phase transition both in the the- 
ory of superfluidity and superconductivity and in economics (a stock price 
break-down, a default, etc.), and, quite unexpectedly, it turned out that, in 
view of the natural axiomatics (see [H]), the zeroth-order phase transition 
is related to quantum statistics and thermodynamics. This phase transition 
has not been noticed by physicists, and, perhaps, it contradicts their ideas 
that the free energy can be determined up to a constant. 

2.5 Metastable state 

We consider a simple example of semiclassical approximation of the one- 
dimensional Schrodinger equation 

-h 2 y n + u(x)y n = X n y n , y n (x) 6 L 2 , h < 1, 

where u(x) = [x 1 — l) 2 + qx and q > is a constant. 
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In classical mechanics, this equation gives the picture of "two cups with 
a barrier between them." If a particle is at the bottom of the right-hand cup 
with higher walls, then it can get into the other cup with lower walls only if 
the barrier disappears. 

From the viewpoint of the Gibbs postulate, as 6 — > 0, the particle must be 
in the deeper well. Nevertheless, it is obvious that if h <C 1, then the particle 
will stay in the well with higher walls for a very long time. In this case, 
the summation in the Gibbs formula must be performed over the subset of 
eigenvalues corresponding to the well with higher walls whose eigenfunctions 
tend to zero as h — > outside this well. Moreover, the temperature 6 must 
be not too high, so that the eigenvalues above the barrier do not play any 
role in the Gibbs formula. 

Such a state at a local minimum of the potential well u(x) is an example 
of a metastable state. 

If we consider the matrix element of transition from the lower level corre- 
sponding to the well with higher walls to the lowest level Ao (corresponding 
to the bottom of the deep well), then it turns out to be exponentially small 
with respect to h, but any transition to higher levels is forbidden accord- 
ing to the "energetic efficiency" law (i.e., we consider perturbations by the 
"friction"). 

2.6 Superfluidity 

N. N. Bogolyubov used Landau's ideas (see the footnote on p. 219 in [T5] ) to 
show that superfluidity is not a fluid motion, not a dynamics, but a state of 
fluid such as, for example, ice or vapor for water. 

This state corresponds to a metastable state of the system such that 
any transition from this state to the normal state is almost forbidden: it is 
exponentially small as N — > oo. 

N. N. Bogolyubov proved this rigorously under the assumption that the 
system of Schrodinger equations is periodic; in other words, the Schrodinger 
equations were considered on a torus. The spectrum of superfluid velocities 
(the energy levels corresponding to the relevant momenta) was discrete. This 
readily distinguishes the state of superfluidity from the hydrodynamics of 
fluids. In the limit as the torus radius tends to infinity, the spectrum does 
not become, as usual, a continuous spectrum, but becomes an everywhere 
dense point spectrum. 

At his time, the author also made such a mistake studying a problem for 
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the one-dimensional Schrodinger equation. He assumed that the spectrum 
was continuous, while the spectrum turned out to be an everywhere dense 
point spectrum. This concerns the counterexample to the Kostyuchenko- 
Morin-Gelfand-Shilov theorem given by the author in 1960. All of them and 
all the specialists agreed with this counterexample and found mistakes in 
their own works. Only 20 years later, Harris found a "hole" in the author's 
proof, which was announced by Barry Simon in his extensive work. After 
this, S. Molchanov and the author found that the example was faulty and 
the spectrum was not continuous, as it seemed to be obvious to everybody, 
but it was an everywhere dense point spectrum. Only in 1987, the author 
and S. Molchanov, using the same idea, could find a counterexample with a 
purely continuous spectrum [To] . 

In the author's opinion, the fact that the spectrum is everywhere dense 
pointwise in the limit can be easily explained from the physical viewpoint. 

Indeed, if the system is in a state with a superffuid velocity v, then 
its transition to a state with larger velocity is forbidden by the "energy 
state efficiency" law, and its transition to a state with any lesser velocity 
is forbidden by the fact that any decrease in the velocity contradicts the 
notion of superfluidity. From the mathematical viewpoint, this means that 
the larger the torus radius, the less (exponentially less) is the transition from 
one state to another. 

2.7 Zeroth-order phase transition 

1.7.1. The author explained the spouting effect discovered in 1938 by 
J. Allen and H. Jones when a superfluid was "flowing" through a capillary 
of diameter 10 _4 cm (in fact, this superfluid was at the superfluid level of a 
metastable state). 

The author used the two-level model to show that, at a point heated (by 
light) till the phase transition temperature 9 C , the heat capacity becomes 
infinite, the entropy has a jump, and the free energy decreases to its lower 
value, i.e., to the point entered by the curve issuing from the ground state 
heated to the temperature 8 C . This means that a zeroth-order phase transi- 
tion occurs. In the present paper, we show that the same picture also appears 
in N. N. Bogolyubov's model of a weakly nonideal Bose gas [T7j . 

This phenomenon can be easily explained if we assume that superfluidity 
is not a thermodynamic state, but the motion of a fluid. At 9 = 9 C , the fluid 
becomes viscous and cannot penetrate through a thin capillary. 
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But the point is that this is not any motion but a state, and then this is 
a zeroth-order phase transition. 

The following question arises. What will happen to superfluidity not in 
a capillary, but in a rather thick pipe? Where is the zeroth-order phase 
transition? At 9 = 9 C) the fluid passes from the superfluid state to the state 
of the usual fluid with viscosity and begins to flow according to usual laws of 
hydrodynamics. After a while, the motion stops, and the transition from the 
superfluid state to the state of a fluid at rest is a thermodynamic transition. 
The total intermediate flow is the hydrodynamics, which must be neglected 
on our time scale. 

1.7.2. Let us consider the example studied in Sec. 1.5 in more detail. 
We slowly vary the constant q. We note that, in thermodynamics with the 
field taken into account, there are two more thermodynamic quantities: the 
field strength and the charge. Thus a variation in q is a variation in a 
thermodynamic variable. 

We show that the passage of q through the zero point results in a zeroth- 
order phase transition. Indeed, the resonance occurs at q = 0: the eigenfunc- 
tions are already not concentrated in one of the wells and the probability 
(the square of each eigenfunction) is identically distributed over both wells. 
The number of eigenvalues is "doubled" and the Gibbs distribution has a 
jump. Precisely here, one can see the role of time. Some time is required for 
half the function, decreasing outside the first well, to be "pumped" into the 
second well. But the Gibbs formula does not take this into account. It may 
happen that we shall wait for this transition for a long time, as it was in the 
preceding example with a thick pipe. 

If q becomes negative, then all the eigenfunctions remain both in the 
first and in the second well, and the Gibbs formula is taken over all the 
eigenvalues, rather than over some of them. 

1.7.3. Finally, we consider the effect of transition into the turbulent flow 
for fluid helium, which, in fact, is very close to Landau's idea concerning the 
energy pumping between large and small vortices. 

As will be shown in another paper, just the resonance that occurs be- 
tween vortices of these two types results in a zeroth-order phase transition, 
which sharply changes the thermodynamic parameters from the viewpoint 
of thermodynamics in which we take into account the large scales of time 
between transitions and, naturally, the averaging over these times. 

1.7.4. Since only the values of bank notes are important, while their 
numbers do not play any role, and the interchange of two bank notes of the 
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same denomination is not an operation, the bank notes obey the Bose statis- 
tics. The averaging of gains is a nonlinear operation, as well as the addition. 
As was shown above, the only nonlinear "arithmetics" (the semiring) that 
satisfies the natural axioms of averaging for Bose particles (bank notes) has 
the form a © h = ln(e a + e b ). 

This lead to a formula of Gibbs type. A variation in (3 = 1/T can be 
treated as a variation in the rouble purchasing power caused, for example, 
by printing a lot of new bank notes. After this, for a period of time, the 
balance (equilibrium) is again established. 

The usual financial efficiency plays the role of energetic efficiency. The 
zeroth-order phase transition is either a default or a crisis |14[ I18j. 



3 An exactly solvable model 

First, we consider the one-dimensional Schrodinger equation for a single par- 
ticle on the circle 

Hip k (x) = Etjj k (x), ip k (x - L) = i/j k (x), (1) 

where ipk( x ) is the wave function, x takes values on the circle, and if is a 
differential operator of the form 

H = e(-ihd/dx), e(z) e C°°. (2) 

Here h is the Planck constant, ip k { x ) — ^xp(ip k x), where p k = 2nhk/L, k is 
an arbitrary integer, and the corresponding eigenvalues are equal to 

E k = e(p fc ). (3) 

We pass from the Hamiltonian function e(p) to a discrete function 7{p) of the 
form 

e(p) = e(nAp) for Ap(^n — - j < p < Ap^n + -J , (4) 

where n is an arbitrary integer and Ap is a positive constant. Then the 
Hamiltonian (2) changes appropriately, and we denote the new Hamiltonian 

by H. The eigenfunctions of this operator coincide with ipk(x), and the 
eigenvalues (3) become 

E k = e( Pk ). (5) 
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In what follows, we assume that the constant Ap takes the form 

2nhG 

Ap = — , (6) 

where G is a positive integer. In this case, it follows from (4) that the set of 
energy levels (5) is the set of G-fold degenerate energy levels 

\ n = E Gn = e(p Gn ). (7) 

The Schrodinger equation for N noninteracting particles has the form 

HnVfa, ...,x N )= E^(x u ...,x N ), (8) 

where ^(x±, . . . , x^) is a symmetric function of the variables x±,...,xn 
(bosons). The Hamiltonian H N is given by the formula 

H N = Y^Hj, (9) 

where Hj is the Hamiltonian of the particle with number j, which has the 
form 

t, = «(-*£). do) 

The complete orthonormal system of symmetric eigenfunctions of the Hamil- 
tonian (9) has the form 

* { jv}(zi, ...,x N ) = Symm XuXN ip {N} (xi, ...,x N ), (11) 

where Symm N is the symmetrization operator over the variables x±, . . . , xn, 
{N} is the set of nonnegative integers N k , k e Z, satisfying the condition 



N * = N > ( 12 ) 



fe=— oo 



and the function ip{ N y(xi, . . . , xn) is equal to 
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ip {N }(xi,...,x N ) = Y[^k s (x s ). (13) 
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Here the indices k±, . . . , are expressed in terms of the set {N} using the 
conditions 



N 

k s < k s+1 for all 1 < s < N - 1, 5 kks = N k for all k G Z, (14) 

and S k k' is the Kronecker symbol. The eigenvalues of the Hamiltonian (9) 

are 

CO 

E({N}) = E k N k . (15) 

k=— oo 

We consider the interparticle interaction of the following form. We assume 
that the particles interact in pairs and the interaction operator for particles 
with numbers j and k has the form 

V jk = -^W{Hj-H k ), (16) 

where V > is the interaction parameter and the function W(l;) is given by 
the formula 

f 1 for If I < D, 
W ® = \0 Jl&>D. (17) 

Here D > is the parameter of the width of interaction with respect to en- 
ergy. The operator (17) corresponds to the interaction under which particles 
in a pair attract each other and radiate the energy quantum —V/N if the 
difference between their energies is less than D and do not interact at all if 
the difference between their energies is larger than D. The Hamiltonian of 
the system of N particles with interaction (17) has the form 

N ^ N N 

j=i j=i k =j+i 

In view of (16), the sums in (18) commute, and hence the set of eigenfunctions 
of the Hamiltonian (18) coincides with (11). It also follows from (16) that 
the corresponding eigenvalues are equal to 

OO t , oo oo 

E(W)= E EkN k -— W(Ek — Ei)(N k Ni — 5 k iN k ). (19) 

k=— oo k=— oo l=— oo 
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In what follows, we assume that the interaction width is sufficiently small 
and satisfies the condition 



D < min |A n - A m |. (20) 

The set of energy levels E k) k e Z, coincides with the G-fold degenerate set 
of levels (7). Hence, by (20), the set of energy levels (19) of the system of N 
particles under study can be written as 

oo . , oo 

E({iV}) = E - E Kn(N n - 1), (21) 

n=— oo n=— oo 

where the level E{iV}) has the multiplicity 

r( M) = (G + ^-. 1 )' . (22) 

Here {N} denotes the set of nonnegative integers JV„, n G Z, satisfying the 
condition 

oo 

E N n = N. (23) 

n=— oo 

Let us consider the statistical sum for the system of N bosons with Hamil- 
tonian (18). Since the energy levels and their multiplicities are given by for- 
mulas (21) and (22), respectively, the statistical sum at temperature 9 takes 
the form 

Z(9,N) = Er({iV})exp((-E{iV})^). (24) 

{N} 

Here the summation is performed over all sets {N} with condition (23) taken 
into account. 

In what follows, we assume that G depends on N and the following con- 
dition is satisfied: 

G 

lim — = g > 0. (25) 

N ^oo N 

By F({N}, 9) we denote a function of the form 

F({N},9) = E({iV}) -01n(r({jV})), (26) 
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and by {N } the set of nonnegative numbers N®, n e Z, for which the 
function (26) is minimal under condition (23). 

Now we consider the problem of finding the minimal value of the func- 
tion (26) under condition (23). In the limit as N — > oo and under condi- 
tion (25), the point of minimum has the form 

N n (Q) — N(m n (9) + o(l)), (27) 
where m n (9), n G Z, is determined by the system of equations 

\ n -Vm n + 9ln( mn )=u(0), neZ, (28) 
\g + m n J 

and fi(8) can be found from the equation 

oo 



m « = L ( 29 ) 



n=— oo 



The substitution of (27) into (26) and then the use of the asymptotic Stirling 
formula give the following relation for the specific free energy: 

f{e) = £m f(9,N)=hm ^ 

N~ >oo N^oo iv 



^2 (^nfn n - + ^6>m n In J -9{g + m n ) In ^1 + 
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(30) 



where, for brevity, we omit the argument 9 of m n (9). We introduce the 
notation 

OJn = A n - ^m n . (31) 
In the notation (31), the system of Eqs. (28) and (29) takes the form 

u n (9) = \ n - V - 9 neZ, (32) 

exp((uj n - n)/9) - 1 

oo 

^ exp(K-/i)/#)-l = L (33) 

n=— oo vv 

The system of Eqs. (32) and (33) exactly coincides with the temperature 
Hartree equations ^H] for the system of N bosons with Hamiltonian (18). 
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We shall study the solutions of Eqs. (28) and (29). For = 0, the system 
has many solutions which we number by the integer /: 

m<? = 5 ln , n,leZ. (34) 

Among all the numbers /, we choose those that satisfy the condition 

v nl = X n - \ + V > for all n ^ I. (35) 

For these numbers, there exist solutions of Eqs. (28) and (29) converging 
to (34) as 9 — > 0. The asymptotics of these solutions as 9 — > has the form 

Vn^Z, 1-mf^^exp . (36) 

Thus, depending on the spectrum A n , n G Z, for sufficiently small values 
of the temperature 0, the system of Eqs. (28) and (29) has many solutions. 
These solutions, along with the point of global minimum, also contain points 
of local minimum of the function (26). The values of the function (26) at the 
points of local minimum are equal to the free energy of metastable states. 
We consider the function (26) at 9 = 0. In this case, it coincides with the 
energy spectrum of system (21). We consider the energy of the system for 
the case in which almost all particles are at the energy level A/, which means 
that following conditions holds: 

N n < N Vn^l. (37) 

Deriving N t from relation (23) and substituting the result into (21), we see 
that, in view of (37), the energy spectrum of the system in this domain has 
the form 

E(iV) » X t N - ™ + ^(A n - X t + V)N n . (38) 

To the Hamiltonian (18), there correspond the Hartree equation and the 
system of variational equations. To each I e Z, there corresponds a solution 
of the Hartree equation, and this solution describes the state 

N^=N5 nl . (39) 

Moreover, the eigenvalues of the system of variational equations correspond- 
ing to this solution of the Hartree equation coincide with 

Vni = A„ - Xi + V, n^l. (40) 



m 



(i) 



~ 9 ex P - 



Vnl 
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In pm I2T] , it was shown that if the eigenvalues of the system of variational 
equations for the solution of the Hartree equation are real and nonnegative, 
then such a solution corresponds to the ground state or a metastable state 
of the system. This means that to all I for which condition (35) holds at 
9 = 0, there correspond metastable states of the system of bosons. As in the 
case of zero temperature, to solutions of Eqs. (28) and (29) for 9^0, there 
corresponds a temperature metastable state if the point (27) is a point of 
local minimum. Now we note that, for very large temperatures, the system of 
Eqs. (28) and (29) has only one solution corresponding to the global minimum 
of the function (26). The asymptotics of this solution as 8 — > oo has the form 

—xjo' (41) 

The uniqueness of the solution at large temperatures means that all metastable 
states disappear with increasing temperature. The temperature at which a 
metastable state disappears is critical for this state. 

Let us consider the behavior of the entropy and the heat capacity of 
metastable states when we approach the critical temperature. We consider 
the metastable state to which there corresponds a solution of Eqs. (28) 
and (29) converging to (34) as 9 — * for some I for which (35) holds. In what 
follows, we assume that the function e{p) satisfies the condition e(0) < e(p) 
for p ^ 0. Then, according to (7), we have A < A; for I ^ 0. Hence the 
solution of Eqs. (28) and (29) converging to (34) as 9 — > for I = corre- 
sponds to the temperature ground state of the system of N bosons. Moreover, 
condition (35) becomes equivalent to the condition A; — A < V. 

We assume that this inequality holds for I ^ 0. The condition that the 
corresponding solution m^iO) of Eqs. (28) and (29) determines the point of 
local minimum of the function (26) under condition (23) by formula (27) can 
be written as the system of inequalities 

ai " ws - y+ ,,i;> w( ;: m .o (())) >° <«> 

°rv)=-v + „, „, <o, -e4^<l 



We note that these inequalities hold for (36) as 9 — ► 0. Inequalities (42) 
follow from the condition that the second variation of the function (26) is 
positive; the variation is calculated under condition (23). 
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The metastable state disappears at the temperature at which the last 
inequality in (42) becomes an equality. We denote this critical temperature 
by #. It follows from (42) and Eqs. (28) and (29) for 9 < # that mi\9) is 
an increasing function of the variable 9 for n ^ I and mf\9) is a decreasing 
function of 9. Moreover, we see that mf ] (9) > m { n ] {9) > m n , (9) if A n < X n i 
and n, n' ^ I. 

From (30) we obtain the following expression for the specific entropy of 
a metastable state in the limit as N —>■ oo: 

* W W= E ((^ + m«)ln(l + ^ )- m (0l n (^L , (43) 

n=-oo ^ ^ 9 / \ 9 / J 

where, for brevity, we omit the argument 9 of rrin (9). Differentiating (43), 
we obtain 

ds ftng) J g + m ff m (0 \ 

The last inequality follows from the properties of rrin\9). Since the last 
inequality in (42) becomes an equality at the critical temperature, we can 
show that, as 9 — > 9c — 0, the solutions of Eqs. (28) and (29) behave as 
follows: 

m^{9) - m«\9P) « ^Sy- Vn e Z, (45) 

Oi n \9c ) 

where is a negative number. The substitution of (45) into (43) shows 
that the derivative of the specific entropy with respect to the temperature 
(this quantity is equal to the heat capacity divided by the temperature) tends 



to infinity as 9 — ► 9^ — according to the law 1/y 9? — 9. This means that 
the projection of a Lagrangian manifold, corresponding to a metastable state, 
on the #-axis is not uniquely determined in a neighborhood of the critical 
temperature 122] • The derivative of the temperature with respect to the 
entropy vanishes as the critical temperature is approached. Therefore, as well 
as in view of (44), as was already pointed out, the Lagrangian manifold is 
uniquely projected on the s-axis. We note that it follows from the properties 
of rrin\9) that the following conditions hold for 9 < 9c^: 

mf{9)<m${9) Vn^O,/, (46) 
o${9) < a n l \9) Vn=£0,l. 
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Inequalities (46) mean that the potential barrier between the energy levels A/ 
and Ao is less than the energy barrier between A; and A n for n ^ 0,1. This 
means that, as the critical temperature is attained, the system of bosons 
under study which is in the metastable state with number I changes its state 
in a jump and passes to the temperature ground state. This is a zeroth-order 
phase transition, since not only the entropy and heat capacity, but also the 
free energy have jumps. 

If the quantity min n ^ |A„ — Ao| = 5 is small, then the difference be- 
tween the free energies in the zeroth-order phase transition from the lowest 
metastable state to the ground state is also small, and the heat capacity in 
this transition has a singularity. The asymptotics of the statistical sum near 
the critical point is given by the asymptotics of the canonical operator in 
a neighborhood of a focal point. This asymptotics has the form of an Airy 
type function of an imaginary argument and has a singularity as iV —>■ oo. 
The parameters i5 C 1, JV > 1, and L ^> 1 can be chosen so that the form 
of the phase transition point coincides exactly with the A-point |23*]-[33j. 

4 Critical temperature and the phase transi- 
tion for a weakly nonideal Bose gas 

Following the Bogolyubov paper [T7j, we consider the system of N bosons on 
a three-dimensional torus with sides of length L. Suppose that the bosons 
interact pairwise and eV(x — y) is the potential of this interaction, i.e., a 
function on this torus, and e—l/N. At the temperature 9, the Bogolyubov 
variational method gives the following expression for the free energy of such 
a system 1 : 

F(e,k ) = F(9, feo, {N(6)}), (47) 

1 Thc same equations arc obtained from the asymptotics of ultra-secondary quantized 
equations of thermodynamics |19|. 
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where F(9, {N}) is the function 

h 2 k 2 e ~ ^ / h 2 k 2 h 2 h 2 \ 

k^ko v 



+ 213 E (^(0) + V(*-A/))JV fc i^ 

k^k'^ko 



(48) 



+ ^ E ^(°)^ + 2ls^(°) ((! - 2iV ) E ^ 
+ £ N k N k , + 2Y,Ni) 

k^k'^ko k^k ' 

+ 6j2i N k HN k ) - (N k + 1) \n(N k + 1)) . 

k=ik 

Here m is the mass of bosons, H is the Planck constant, k and fc' are arbitrary 
three-dimensional vectors of the form 27r/ L(ni, n,2, n^), ni,ri2, n 3 are integers, 
k is a distinguished vector of this form, and V is given by the formula 

Finally, {N} denotes a set of variables N k) k ^ k , and {N(9)} is the mini- 
mum point of the function (48). We note that the function (48) corresponds 
to the case in which the system of bosons contains a condensate consisting of 
particles with momentum hk . Suppose that k ^ is less than the Landau 
velocity, i.e., ko must satisfy the condition 

fh 2 k 2 h 2 B V(k-ko)\ 

min + V , °M > 0. 50 

k^k \ 2m 2m L 3 J 

To such a ko, there corresponds a metastable state. The equations for deter- 
mining the minimum point of the function (48), 
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are equivalent to the system of the temperature Hartree equations. In view 
of (50), Eqs. (51) have solutions for sufficiently small temperatures 9, and 
the asymptotics as 9 — > of these solutions has the form 

N k « e -"W e VA; ^ fco, (52) 

where _ 

ft 2 * 2 h 2 k$ v(k - fc ) 
= - + ' (53) 

We note that the substitution of (52) into (48) leads to an expression different 
from the free energy obtained by the summation over the quasiparticles. 
Suppose now that ko is less than the Landau velocity, but still is very close 
to it. This means that the minimum in (50) is attained at k — 0, and this 
minimum tends to as N — > oo. To be definite, we assume that 

^ofco = ^r, (54) 

where 1/2 < a. Since Ukk for /c 7^ is a quantity of the order of 1, it follows 
from (54) that 

N k < N (55) 

at any temperature. In view of (55), in the limit as N — > 00, the leading 
contribution to the function (48) is given by the terms containing N ; thus, 
the terms that do not contain N can be neglected. In this case, the system 
of Eqs. (51) for determining the minimum point becomes the equation 

where 

P = 2(2V(k ) + V(0)) > 0. (57) 

Equation (56) can trivially be reduced to the quadratic equation. An analysis 
of this equation shows that it has two solutions for < 9 < 9 C , where 

B ^im^' (58) 

and Eq. (56) does not have solutions for 9 C < 9. 



21 



As 6 — > 6 C — 0, the solutions of Eq. (56) have the form 



The substitution of (59) into (48) gives the free energy whose second-order 
derivative becomes infinite as 6 — > 9 C — 0. This means that the heat capacity 
of a metastable state increases significantly near the temperature that is 
critical for this metastable state: 

cvm = -^l? o) » , 1 ■ (60) 

Thus, it turns out that the properties of metastable states in the Bogolyubov 
model are similar to the properties of metastable states in the two-level 
model. Namely, metastable states disappear at the critical temperature, and 
the heat capacity of these states increases according to the law 1 / \/8 c — 9 as 
the critical temperature is approached. 

Hence we have the same picture of the zeroth-order phase transition and 
the spouting effect for a weakly nonideal Bose gas in the Bogolybov model. 
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